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Abstract — The first part of this paper presents a simple and 
systematic tec hni que for constructing multidimensional A/ -ary 
phase shift keying (MPSK) trellis coded modulation (TCM) codes. 
The construction is based on a multilevel concatenation approach, 
in which binary convolutional codes with good free branch 
distances are used as the outer codes and block MPSK modu- 
lation codes are used as the inner codes (or the signal spaces). 
Conditions on phase invariance of these codes are derived and 
a multistage decoding scheme for these codes is proposed. The 
proposed technique can be used to construct good codes for both 
the additive white Gaussian noise (AWGN) and fading channels 
as is shown in the second part of this paper. 

Index Terms — Concatenation, MPSK modulation, multidimen- 
sional trellis coded modulation, multistage decoding. 


I. Introduction 

S INCE the publication of the celebrated paper by Unger- 
boeck on trellis coded modulation (TCM) [1], there has 
been a boom of research in this area. Over the last fourteen 
years, researchers have proposed various techniques of con- 
structing modulation codes using both convolutional codes 
(TCM) [l]-[7] and block codes [block coded modulation 
(BCM)] [8]-[14]. Almost all existing techniques for con- 
structing TCM codes rely heavily on computer searches to 
find good TCM codes. These techniques work very well for 
small code complexities and rates. However, for large code 
complexities and high rates, the search becomes extremely 
time consuming (if not impossible) and a more systematic 
technique of construction is required. Most of the problems 
associated with the algebraic construction of TCM codes arise 
due to the lack of in-depth knowledge of convolutional codes. 
In addition, the nonlinearity of the mapping function (true for 
most signal constellations) which maps the coded output bits 
of the convolutional encoder onto the signal set, complicates 
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the problem further. BCM codes on the other hand, have 
the advantage of being extremely rich in algebraic structure 
and phase symmetry, as has been shown in [10]— [ 13]. BCM 
codes however, have the disadvantage of being slightly poor in 
performance for low signal-to-noise ratio (SNR), as compared 
to TCM codes of the same decoding complexity, due to the 
large number of nearest neighbors. 

Pietrobon et ai extended Ungerboeck’s results to multidi- 
mensional MPSK signal constellations [3]. They proposed a 
set partitioning technique for multidimensional MPSK sig- 
nal constellations similar to Ungerboeck’s set partitioning 
technique and then used computer search to design multidi- 
mensional MPSK TCM codes. However, due to the limitations 
of computer search, as were outlined above, they restricted 
themselves to 4 x 2 dimensions. In addition, to reduce the 
search complexity, they placed some other restrictions on 
the computer search. Multidimensional MPSK TCM codes 
have various advantages over two-dimensional (2-D) Unger- 
boeck TCM codes, the main ones being: 1) higher spectral 
efficiencies can be achieved, 2) codes constructed over mul- 
tidimensional MPSK signal constellations have better phase 
invariance properties than that of 2-D Ungerboeck MPSK 
codes, and 3) lower average decoding complexities to achieve 
the same performance. 

A common point to be noted among all the construction 
techniques available in literature (whether TCM or BCM) is 
that the modulation codes constructed by these techniques 
require large decoding complexity to achieve large coding 
gains. The large decoding complexity of these codes makes 
them impractical for applications where high reliability and 
high data rates are required. As such, what is required is a 
multistage decoding technique which reduces the decoding 
complexity, while maintaining good performance. 

This paper presents a simple and systematic technique for 
designing multidimensional MPSK TCM codes with minimal 
computer search. The technique will be used to construct good 
codes for both the AWGN and fading channels. Though the 
main emphasis has been to construct codes for the MPSK 
signal constellation, the results are applicable to other signal 
constellations as well and modifying the existing construction 
for other signal constellations is straight forward. This paper 
is organized as follows: Section II of the paper presents a new 
concept, branch distance of convolutional codes, which will 
be used extensively in the later sections. Section III outlines 
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the basic construction technique of the proposed codes, and, 
in addition, shows that the codes constructed in [ 3 ] turn out 
to be a special case of the proposed construction. Section 
IV discusses phase invariance. In Section V, a multistage 
decoding algorithm for the proposed codes is presented and 
it’s decoding complexity is discussed. Section VI concludes 
by discussing the design rules for constructing good codes 
using the proposed technique. 

n. Branch Distance of Convolutional Codes 

For two code sequences « and v in a binary linear con- 
volutional code, the branch distance between them, denoted 
is defined as the number of branches in which t* and 
v differ (or equivalently, this is simply equal to the number of 
nonzero branches in t*0 v, where © denotes binary addition). 
For a code sequence u in a binary linear convolutional code, 
the branch weight of u denoted lUfc(ti) is simply the number 
of nonzero branches in u (or equivalently wt,(u) is the branch 
distance between u and 0, where 0 refers to the all-zero code 
sequence, i.e., u/&(ti) = dfe(ti, 0 )). The minimum free branch 
distance of a convolutional code C, denoted d B . free, is the 
minimum branch distance between any two code sequences, 
i.e., 

^B-free = min{cA(ii, v): u,v e C and u^v). (2.1) 

Theorem 1 : For a rate k/n feedforward binary linear con- 
volutional code of total encoder memory 7, its minimum free 
branch distance, ds-free> is upper bounded by 1 + [y/k\. 

Proof: Let the k inputs to the encoder be denoted as 
A, A, • ■ * , h and let the encoder memories associated with • 
input A be 7^ for 1 < z < k. Let min{u;fc(ti)} denote the min- 
imum branch weight among all the code sequences associated 
with the binary linear convolutional code. Let min^Lj 7* = 
7 j. Consider that the binary sequence ( 1 , 0 , 0 , * • •) is fed 
into the input I 3 and the all zero sequence (0, 0 , 0 , ••■) is 
fed into the remaining inputs. The branch weight of the 
resulting code sequence is upper bounded by 1 4- 7 j. Hence, 
minjtt^tt)} < 1 + 7^. Since the code is linear, this also 
corresponds to an upper bound on the minimum free branch 
distance, i.e., d B -free < {1 + min*^ 7,}. Given any 7 and 
k , the idea is to maximize ds- free- Hence, max^fc^B-free) < 
max 7> ]fc{l -t-minfLj 7 *}, i.e., the best d B -f re€ for a given 7 and 
k is < {1 -f max 7i fc{min*L 1 7*}}. It is readily seen that the 
value of max 7 ,*{min* =1 7*} is [7/fcJ. A A 

Theorem 2 : If d B -free = 1 + [l/k\, then iV B -f re e, the 
number of codewords with branch weight d B - free, is lower 
bounded by (2 P - 1) where p is the number of inputs of 
the convolutional encoder which have an encoder memory of 
[7/fcJ associated with it. 

Proof: Let e x denote the binary sequence ( 1 , 0 , 0 ,***) 
i.e., 1 followed by the all zero sequence and let eo denote the 
all zero binary sequence ( 0 , 0 , 0 ,***). Consider any nonzero 
code sequence u. Then > 1 -f [7/fcJ. Let the p inputs 

which have an encoder memory of [7/fcJ associated with it 
be Ij for 1 < j < p. Consider that eo is fed into the inputs 
Ij for p + 1 < j < k. Also, consider that the inputs Ij 


TABLE I 

Optimum Branch Distance Rate 1/2 Codes 



t Total encoder memory. 

* Minimum free branch distance. 

A Number of codewords with branch distance <f B -free* 
u Free Hamming distance. 

* Number of codewords with Hamming distance ^H-free 
Note: The code generators have been listed in octal, where the 
octal representation of xyz is 4 • x 4* 2 y + z and x, y 
and z denote three binary bits. 


for 1 < j < p can take only one of the two sequences eo 
or e\. Then the convolutional encoder under this constraint 
has (? p - 1) distinct nonzero input sequences. Each of the 
(2 P - 1) sequences will have branch weight < 1 + [7/fcJ. 
Since ds- free = 14- [l/k\ , each of the ( 2 P — 1) sequences thus 
has branch weight 1 + ll/k\ . Hence, N&. f ree > (2 P - 1). A A 

A binary linear feed-forward convolutional code is said to 
be optimal in terms of branch distance if it achieves the upper 
bound as stated in Theorem 1 for a given 7 and k. Also, 
a code is said to be optimal in terms of the free Hamming 
distance, dH-free> if it achieves the maximum free possible 
for a given 7, k and n as specified in [ 15 ]. Note, from Theorem 
1, for a given d B -free, higher encoder memory is required to 
achieve the same dB-free as k increases, i.e., given a certain 
fixed d B -fre« ? there is a tradeoff between complexity and rate. 
In addition, as is shown in Theorem 2 , N^-free also increases 
as the rate increases and hence there is also a tradeoff between 
rate and performance. A search has been performed on rate- 
1/2, - 2/3 and - 3/4 codes to find the best ones in terms of d B -f re * 
and A^-free* The results are given in Tables I-III. 

An important point to note is that codes optimum in terms 
of branch distance may not be optimum in terms of the free 
Hamming distance dn.free and vice-versa. For small values 
of 7, it has been observed that codes optimum in terms of 
branch distance are also optimum in terms of dH-free> however, 
the same does not hold for higher values of 7. From Table 
I, we notice that up to 7 = 7 , the search yields codes 
which meet the upper bound in terms of d B -free> however 
from that point on, the best codes start falling short of the 
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TABLE II 

Optimum Branch Distance Rate 2/3 Codes 
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t Total encoder memory. 

♦ Minimum free branch distance. 

A Number of codewords with branch distance dg.f ree . 
u Free Hamming distance. 

* Number of codewords with Hamming distance ^H-free 
Note: The code generators have been listed in octal, where the 
octal representation of xyz is 4 • x + 2 - V 4- £ and x, y 
and a denote three binary bits. 


TABLE III 

Optimum Branch Distance Rate 3/4 Codes 
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t Total encoder memory. 

* Minimum free branch distance. 

A Number of codewords with branch distance <%.f ree - 
u Free Hamming distance. 

* Number of codewords with Hamming distance djj_f ree . 
Note: The code generators have been listed in octal, where the 
octal representation of xyz is 4 • x + 2 y + z and x, y 
and c denote three binary bits. 


upper bound by one. Codes shown in Tables II and HI 
meet the upper bound, however as the complexity increases, 
N &- free al so starts increasing. Also listed in the tables is the 
dn- free and TVn-fre*, the number of codewords with cfn-free* 
The code generators in the tables have been listed in octal 
with the lowest degree on the left and the highest on the 
right, e.g., (622) 8 = 1 + D + D 4 + D 7 . As an example, 
consider the eighth code listed in Table I. This is a rate- 1/2 
convolutional code with generators 1 -F D + D 4 + D 7 and 
1 + D 2 + D 3 + D A + D 5 + D 6 + D s and d B -free = 8. 


III. Construction of Multidimensional MPSK Codes 

The proposed multidimensional MPSK codes are con- 
structed using a q level concatenation approach as shown 
in Fig. 1. Outer codes in the multilevel concatenation may be 
either block or convolutional, binary or nonbinary. However, 


in this paper we will focus on binary convolutional codes as 
the outer codes. 

Outer Codes: The outer code, C;, at the zth level for 
1 £ * < q is chosen to be a convolutional code of rate ki/rii 
with optimum branch distance for the given rate and state- 
complexity. The parameters k x and n; depend upon the choice 
of the inner codes, as will be clear after the discussion of inner 
codes. Each outer code is selected from the tables mentioned 
in Section II. The reasons for selecting an optimum branch 
distance convolutional code will be clear when discussing 
Theorems 4-6. 

Inner Codes: Let S denote the two-dimensional MPSK 
signal constellation which consists of 2* signal points. Let S m 
denote the set of all m-tuples over 5, where m is a positive 
integer. Since 5 is a two-dimensional signal space, S m is an 
m x 2-dimensional signal space in which each signal point 
is a sequence of m MPSK signals. To construct the proposed 
codes, the signal space is chosen as a subspace of S m , denoted 
Ao- In this paper, Ao is constructed using the multilevel coding 
method proposed by Imai and Hirakawa [8]. 

Using the set partitioning approach proposed by Ungerboeck 
in [1], each signal point in the set S is labeled by a string 
of symbols from GF(2). Since 5 contains 2 e signal points, 
we shall consider a labeling whose set of label strings is of 
the following form: L = {ai 02 ■ • - a*: a x € GF(2) for 1 < 
i < £}. Let A denote the one-to-one mapping from L to 
S. If 0 x 02 •• -a* is the label for a signal point s, then 
$ = A(aia 2 • • - a*). Define an addition “+” on the label set 
L as follows: For two labels, aia 2 • * * ai and • • • a' t , in 
L, aia ,2 • • • at -f a[a f 2 •••©£ = a” a % • * • a" where a" = a,6a- 
for 1 < i < i and 0 is the modulo-2 addition. With this 
addition, L is simply the vector space of all /-tuples over 
GF(2). We call L the label space for 5. 

For 1 < i < /, let Co,; be a binary (m, koj, 5 0 ,;) linear block 
code of length m, dimension kos and minimum Hamming 
distance <5o,;. Let 

V{ = v <fm ) (3.1) 

be a code word in C 0 ,; for 1 < i < l. We form the following 
sequence: 

Vi * V 2 * • • * * V £ 

a / 

,m ). (3.2) 

For 1 < j < m, we regard vijV 2 j * ■ * vtj as the label for a 
signal point Sj in the MPSK signal set S. Then V P i*V r 2 **V r f 
is simply an m-tuple over the label set L and 

A(Vi *V 2 ---*Vt) 

= (A(vi,i V2,l - • • W/,l), \(Vi.2V 2 ,2 ■ ■ ■ Vl, 2 ), 

* ' * , A(Vi im t72,m * ‘ * ) ) 

= (Sl,S 2 , ( 3 - 3 ) 

is an m-tuple over the MPSK signal set S (a sequence of 
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Inner Code 




Fig. 1. A multilevel concatenated TCM system. 


m MPSK signals) which is a signal point in the m x 2- 
dimensional signal space S m . From codes Co for 1 < i < t, 
we form the following set of m-tuples over the label set L : 

Co.i * ^ 0,2 * * * • * 

= {Vi * V 2 * - - - * V t : V 1 e C 0 a,V 2 € C 0l 2 , 

• v t e C Qt i}. (3.4) 

We will denote Cq,i * Co , 2 * ■ • • * Co,* by Qq- Then, is 
a vector space (or a linear code) over L (a subspace of the 
vector space of all the m tuples over L , denoted L m ). Qq 
has 2*°‘ 1+fc ° 2+ * vectors. Hence, the dimension of Oo Is 

cr 0 = &o,i + ko ,2 H h A:o,£. Recall, that for 1 < ? < <j, n t 

denotes the number of output coded bits of the convolutional 
encoder at the ith stage of encoding. Choose 

4* 712 H“ ' ’ • + Tlq = ko i + ko,2 + 1- k 0l t — <7o (3.5) 

Suppose each m-tuple in Do is mapped into an ra-tuple over 
the MPSK signal set 5 by the mapping A(-). Then, we obtain 


the following subset of signal points in S m : 

A 0 = A(Oo) 

= {A(V! * V 2 • * V e ): V t € C 0 ,i, V 2 € Co, 2 
■V t &C 

The set Ao is a subspace of S m with dimension ao- This sub- 
space Ao is actually a fowic Mevel block MPSK modulation 
code of length m [8]— [14]. 

The performance of Ao over the AWGN channel depends 
upon the minimum squared Euclidean distance and the number 
of nearest neighbors. The minimum squared Euclidean dis- 
tance of Ao can be calculated using results of [12]. On the other 
hand, the performance of Ao over fading channels depends 
upon the minimum symbol distance, product distance, number 
of nearest neighbors and the squared Euclidean distance to 
a lesser extent [17]. The minimum symbol distance of Ao is 
given by [17] 6% = minf =1 So Suppose, A 0 has minimum 
squared Euclidean distance Aq and minimum symbol distance 
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In the following, the subspace Ao of S m will be used as the 
signal space for constructing multidimensional trellis MPSK 
codes. Before presenting the code construction, we need to 
define a subspace of ft 0 for partitioning ft 0 . For 1 < j < t, 
let Cij, C 2 j, • • • , C qj be a sequence of linear subcodes of 
Cqj such that 

C q j C C g -u C . • • C Cij C Coj- (3.6) 

Let kij be the dimension and 8ij be the minimum Hamming 
distance of C itj for 1 < i < q. Then Cij is an (m, kij,8 iyj ) 
code. For 1 < i < q, we form the following linear code over 
the labeling space L : ft* = Ci t \*Ci £*• * The dimension 

of this code is Oi = + ki ,2 + • ■ • + It is clear that 

for 1 < i < q, 

ft* C ft;_ i. (3.7) 

It follows from (3.7) that fti, ft 2 , * * * . ft 9 form a sequence of 

subspaces of fto and 

ft* c ft*_i c.-.cflicno. (3.8) 

For 1 < t < let 

A, i A(fli). (3.9) 

Then, A i is a subspace of S m with dimension dim(Ai) = 

Let the minimum squared Euclidean distance of A; be A? and 
minimum symbol distance be 6 l H . Equations (3.8) and (3.9) 
imply that Ai, A 2 , ■ ■ ■ , A* form a sequence of subspaces of 
Ao' and 

A* C A q -i C ••• C Ai C A 0 . (3.10) 

Suppose the binary codes, Cij with 1 < i < q and 1 < j < t, 
are chosen such that 

m = _ 1 — < 7 (3.1 1) 

It follows from (3.5) and (3.11) that 

<7\ = n 2 + 713 4 * • * • + Tl q 
<J 2 = ^13 -|- • • * + Tlq 

< 7 q — l = Tlq 
(Tq =0 

ft 0 and its subcodes fti,ft2,*">ty/ are used to form a 
sequence of coset codes [7]. Let U\ * U 2 * * * • * Ue be a 
vector in fto but not in fti. Then U\ * U 2 * * * * Ut 4- fix is a 
coset of Six in ft 0 and U\ * U% * • * * Ut is called the coset 
representative . Recall ni = £f =1 (A;o,i — fci,*). Hence, there 
are 2 ni cosets of fti in fto. These 2 ni cosets of fti form a 
partition of Ho- Let ft 0 /fti denote the set of cosets in fto 
modulo fti- fto/fti is called a coset code. Let [fto/fti] denote 
the set of coset representatives of the coset code fto/fti- Hence 
ft 0 /Qi = [fto/fti] + fti. fti can be further partitioned using 
ft 2 • in the same way as is outlined above. Partitioning each 
coset of fti in ft 0 on the basis of ft 2 , we form the coset code 
ft 0 /fti/ft 2 - Let [Qi/Q 2 ] denote the set of coset representatives 
in the partition ft 1 /ft 2 - Hence each coset in the coset code 


fto/fti /^2 can be written in the form [ft 0 /fti] + [fti/ft 2 ]+ft 2 * 
Proceeding in this manner, we form the following sequence of 
coset codes: 

B\ — fto/fti 
B 2 = fto/fti/ft2 

Bq = fto/ftl/^ 2 / * * * /ft<7 • 

For 1 < i < q, each coset in B t -i = fto/fti/ *•* /fti-i 
consists of 2 n * cosets modulo ft*. These coset codes are used 
as the inner codes in the multilevel concatenation in which B\ 
is used at the first level and B q at the qth level. 

Let ljq and w' 0 be two distinct points in ft 0 . If these two 
points are in two distinct cosets of B\ then the squared 
Euclidean distance between s = A(u> 0 ) and $' = A(cjq) is 
at least Aq. If the two points ljq and cj' 0 are in the same coset 
of B\ but distinct cosets of B 2 , then the squared Euclidean 
distance between s and s' is at least Aj. Generalizing in this 
manner, it is easy to see that if the two points u>q and u' 0 have 
identical coset representatives in Bj for 1 < j <i, but distinct 
coset representatives for Bi then s and s' have a squared 
Euclidean distance of at least A^_ x . Hence, B 1 is the least 
powerful and B q is the most powerful coset code in terms of 
Euclidean distance. 

The same arguments as above will also hold if the minimum 
squared Euclidean distance at each stage is replaced by the 
corresponding minimum symbol distance. 

Encoding of m x 2-Dimensional TCM Code : Encoding is 
accomplished in q stages, as shown in Fig. 1, and for 1 < 
i < q, the ith level encoding is accomplished in two steps: 1) 
at any time instant t, a message of bits is encoded based 
on the convolutional outer code Ci into an n*-bit coded block; 
and 2) the 7 ^ -bit code block then selects a coset from the 
coset code Bi = fto/fti/ •** /ftt- 

The output at the zth level encoder is a sequence of cosets 
from Bi. All the possible coset sequences at the ith level 
form a trellis, and each branch in the trellis corresponds to 
a coset in Bi, and this trellis is isomorphic to the trellis of 
Ci. Let Vj denote a code sequence in the convolutional code 
Ci and let fa denote the mapping from the n l coded output 
bits of the convolutional code to the 2 n% cosets. Hence, fa{v q ) 
denotes the sequence of coset representatives at the ith stage 
of encoding, corresponding to the code sequence v t *. Hence, 
any code sequence in the m x 2-dimensional TCM code can 
be written in the form 

A(<£i(t;i) + <£ 2 (^ 2 ) + h <l>q(vq))- (3.12) 

At every time instant t, the encoder puts out m MPSK signals. 

A very interesting and special case of the proposed codes 
occurs when q = 2 and the second level outer code is left 
uncoded, as shown in Fig. 2. This structure is equivalent to 
the structure used for the construction of the multidimensional 
codes in [3]. A computer search was used in [3] to find the 
convolutional code to be used at the first level. The computer 
search selected a convolutional code which optimized the 
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multidimensional code both in terms of Euclidean distance 
and number of nearest neighbors. 

A multidimensional code is said to be linear with respect 
to binary addition, if for any two code sequences in the 
multidimensional code, U = A(<£i(ui) + <£ 2 (^ 2 ) + + 

$ q (* q )) and V = A(^!(wi) + $ 2 ^ 2 ) + • • • + $ g ( vq )) 

V © V = A((<£i(t»i) + $2(112) + • • • + $ q (uq)) 

+ ($i(vi) + <fo( vt) + • • • + 4 >,(vq))) 

is also a code sequence, where tij and «j for 1 < i < q denote 
output code sequences of the convolutional code encoder C{ 
at the zth level. Linearity of the code (in terms of binary 
addition) simplifies the error analysis and in addition leads 
to a simpler encoder and decoder. The linear structure leads 
to the following theorems on the linearity, minimum squared 
Euclidean distance, and minimum symbol distance of the 
proposed codes. 

Theorem 3: A multidimensional code is linear with respect 
to binary addition, if all the mappings fa, for 1 < i < q are 
linear. 

Proof: Recall, that any code sequence in a multidimen- 
sional TCM code can be written in the form \{<pi{vi) + 
<£ 2 (^ 2 ) + ■•■ + 4 >q{ v q)) where Vj for 1 < i < q denotes 
the output code sequence of the convolutional code C{ at the 
ith level. The proof then follows trivially from the definition 
of linearity. A A 

Theorem 4: The minimum free squared Euclidean distance 
of a coset trellis code at the jth level, for 1 < j < q is lower 
bounded by D? j) hee > A 2 _ x • d^} fTee , where dg kee denotes 
the minimum free branch distance of the convolutional code 
at the jth level, C 3 . 

Proof: Consider two distinct code sequences, U = 
A(^i(t»i) + + •• - + <t>q{uq)) and V = \{ 4 >i{v x ) + 

^2(^2) + *** + 4 > q (vq)), where tij and Vj for 1 < i < q 
denotes two output code sequences of the convolutional code 
Ci at the ith level. Assume that = Vj for 1 < i<j and 
Uj ^ Vj. At a particular time instant let \{lj) and A(u/) 
be the corresponding transmitted signal points for U and 
V , respectively, where uj and uj 1 E fio* Since for 

1 < i<j and uj ^ vj, hence u j and c J have identical coset 


representatives in B{ for 1 < i < j and hence the minimum 
squared Euclidean distance between A(u>) and A(u/) is at least 
A Since Cj has minimum free branch distance 
hence the two sequences uj and vj are distinct in at least 

^B-free branches. Therefore, the squared Euclidean distance 
between U and V is at least A^_ x • d%} hee . A A 

Theorem 5: The minimum free squared Euclidean distance 
of the overall TCM code is lower bounded by D ^ > 

Proof: Consider two distinct code sequences U and V . 
Using the same notation as developed in Theorem 4, consider 
that 11 ^ = Vj for 1 < i<j and that uj ^ vj. Then, 
Theorem 4 gives us the minimum squared Euclidean distance 
between the two sequences. Since j is arbitrary, the minimum 
squared Euclidean distance between the two sequences is 
obtained by taking the minimum over all the q levels, i.e., if 
D 2 (U , V) denotes the squared Euclidean distance between the 
two sequences U and V, then D 2 (U, V) > mini^j^jA^j * 

^B-fre©}* Since U and V are any two sequences, the Theorem 
follows. A A 

Theorem 6: The minimum symbol distance of the overall 
TCM code is lower bounded by <5# > mini<j< 9 {<5^ -d#} {Tee }. 

Proof: The proof is similar to that in Theorem 5, with 
the only difference that instead of minimum squared Euclidean 
distance we now consider minimum symbol distance. A A 

IV. Code Properties 

A. Spectral Efficiency 

At each encoding time instant, fci + &2 + b k q bits are 

fed into the encoder (Fig. 1), and the corresponding output 
is m MPSK signals. Hence the spectral efficiency of the 
m x 2-dimensional TCM code is (k\ + fc 2 + • • • + k q )/m 
bits/symbol. 

B. Phase Invariance 

Phase symmetry of a code is important in resolving carrier- 
phase ambiguity and ensuring rapid carrier-phase resynchro- 
nization after temporary loss of synchronization [2]. It is 
desirable for a modulation code to have as many phase sym- 
metries as possible. Recall, that the proposed multidimensional 
modulation codes are constructed using q convolutional codes 
and q - hi basic ^-level block modulation codes (Fig. 1). The 
phase invariance of the proposed codes is a function of both 
the inner codes and the outer codes. If convolutional codes are 
used at all the q levels , the phase invariance of the constructed 
modulation codes would depend upon the structure of the 
convolutional codes used, and for most cases the constructed 
modulation codes would have no phase invariance. A special 
case of the proposed codes occurs when the outer code at the 
< 7 th level is left uncoded (Fig. 2 shows this special case for 
q = 2). Most of the codes constructed using this special case 
do have phase invariance. Kasami et al derived conditions on 
phase invariance of basic f-level block modulation codes [16]. 
A slightly modified form of the conditions proposed in [16] 
will be applicable to the proposed codes. 
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The following theorem gives the conditions for the proposed 
modulation codes to be phase invariant under rotation for this 
special case. 

Theorem 7: Let Ao = A(Co,i * Co , 2 * * * * * C<m) and let 
A q -i = A(C 9 _i t i * C q - 1,2 * where C 0j i and 

C q -i,i for 1 < i < t are binary linear block codes of 
length m. For 1 < h < l, the multidimensional MPSK TCM 
code is phase invariant under 180°/2 f_h phase shifts if the 
multidimensional TCM code is linear with respect to binary 
addition and 

1 eC q . hk and (4.1) 

C 0 ,h‘Co,h+i---C 0 j-iCC q - hj for h<j<i (4.2) 

where 1 denotes the all-one binary sequence of length m, 
and for two-binary m-tuples a — (ai, a 2 , • * * a m ) and 6 = 
{bi,b 2 , - • -bm)^ * b = (ai * &i,a 2 • 6 2 ,--*,a Tn • 6 m ), where 
ai • 6i , for 1 < i < m denotes the logical product of a t and b { . 

Proof Appendix: 

If the outer code at the q\h level is left uncoded, sequences 
of signal points from are valid code sequences. The 

best phase invariance that can be achieved for the overall 
multidimensional code in this case is equal to the phase 
invariance of A 7 _i. The conditions as stated in Theorem 7 
provide a set of conditions which guarantee a certain phase 
invariance for the overall multidimensional MPSK TCM code 
independent of the convolutional codes chosen. Most codes 
designed using the proposed technique, do achieve the best 
possible phase invariance (i.e., of A^-i). 

V. Multistage Decoding Algorithm 

One obvious way of decoding a TCM code proposed in 
Section III, is to form a super trellis for the code, which is 
obtained by taking the direct product of the trellises of the 
convolutional codes at the q levels. The complexity associated 
with this technique (for most cases) would be tremendous. 
We will focus on a multistage decoding scheme, in which 
the decoding is carried out in q stages, corresponding to 
the q levels of the multidimensional TCM code. Let V = 
(si, s 2 , 53 , ■ • ■) be the transmitted code sequence, where $i 
for 1 < i < oo denotes a signal point in the MPSK 
signal constellation and let R = (i%r 2 ,i% • • •) denote the 
corresponding received sequence. Using (3.12), V can be 

written in the form V = A(<£i(i/i) + ^ 2 (^ 2 ) H 1- <t>q( v q)) 

where vj for 1 < i < q denotes a code sequence in the 
convolutional code Ci. 

First Stage of Decoding : At the first stage, t>i is estimated 
using the received sequence R. Recall, that at the first stage 
of encoding, the trellis is isomorphic to the trellis of the 
convolutional code C\ used at the first level, with each branch 
of the trellis corresponding to a coset in B\. Each coset 
in B\ can be written in the general form ljq 4- fli, where 
wo *= [fio/Oil- Let us ca ^ isomorphic trellis C\. Hence, 
each branch of C\ consists of 2° l points, corresponding to 
the 2 ai points in Dj. The trellis Ci is used to form the trellis 
A(Ci), where 


The trellis A(Ci) will be used for decoding at the first stage. 
Any code sequence in A(Ci) can be written in the form 

\(<f>i(ui) + a/i ) (5.2) 

where ui is a code sequence in C\ and wi is a sequence 
of points from Hi, i.e., u)\ = wi f 3 , • • •): u>\,i G 

f2i for 1 < i < 00 }. Standard soft-decision Viterbi decoding 1 
is performed on R using the trellis A(C 1 ). This yields a code 
sequence X((f>x(vi) +«i) in A(C'i) which is closest to the 
received sequence R in terms of squared Euclidean distance. 
The code sequence t>i forms an estimate of the sequence 
v\. <j)\ denotes a sequence of points from Q\. Since t>! is a 
code sequence in Ci, the estimate of the information sequence 
associated with the first level can be obtained from t>i. 

The ith Stage of Decoding: The second and subsequent 
stages of decoding are very similar to the first stage of 
decoding. For 2 < i < q, let us consider the ith stage of 
decoding. The previous i - 1 stages of decoding give us 
estimates of denoted by i)j for 1 < j < (i — 1). Using 
arguments similar to that given above, we form the isomorphic 
trellis Ci , where any code sequence in Ci can be expressed 
in the general form 

^i(vi) -F <£ 2 (^ 2 ) 4- h <Pi-i(vi-x) 4- <M u i) + u3x 

where Uj is a code sequence in the convolutional code at the 
ith level, Ci and u>i is sequence of points from D*. Each branch 
of Ci consists of 2°' points, corresponding to the number of 
points in Qi. The trellis Ci is used to form the trellis A (Ci), 
where 

A(Ci) £ {A(»):t/eCi}. (5.4) 

The trellis A (Ci) will be used for decoding at the ith stage. 
Standard soft-decision Viterbi decoding is performed on R 
using the trellis A(Ci). This yields a code sequence 

A(<£i(t>i) + <£ 2 (^ 2 ) + F 4- <t>i{vi) + u>i) (5.5) 

in A (Ci) which is closest to the received sequence R in terms 
of squared Euclidean distance, where t>i is a code sequence in 
the convolutional code used at the ith level, C{ , and tdi is a 
sequence of points from Hi. The code sequence t ); forms an 
estimate of the sequence Vi. Since Vi is a code sequence in 
Ci , the information sequence associated with the ith level can 
be obtained from v*. 

The branch metric (squared Euclidean distance) for each 
branch in A(C,),1 < i < q, is calculated by taking the 
m received signals corresponding to that branch and finding 
the element in the coset corresponding to that branch, which 
is closest to the m received signals in terms of Euclidean 
distance. This process of finding the closest element in the 
coset is termed as closest coset decoding . The Euclidean 
distance corresponding to the closest element in the coset 
becomes the branch metric. If m is small, calculation of the 

1 We will use minimum squared Euclidean distance as the decoding metric 
for both the AWGN and fading channels. 


A(Ci) = {A(*):®eCi}. 


(5.1) 
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branch metric does not represent a formidable task, however 
if m is large and if ft;, 1 < i < q, has trellis structure then a 
trellis can be used to calculate the branch metric. In addition, if 
the number of states associated with the trellis structure of ft; 
is big, multistage decoding for ft; can be used to further reduce 
the decoding complexity. Multistage decoding of ft* would be 
carried out in the same way as proposed in [ 10 ] and [ 11 ]. 

Another way of reducing the decoding complexity associ- 
ated with closest coset decoding would be as follows: Consider 
a trellis C* up , where any code sequence in the trellis C* up can 
be written in the following form: 

4>\{vi) + < t > 2(62) + • * * + + M u i) + w i UP ( 5 - 6 > 

where is a sequence of points from J?® up , and the rest of 
the sequences are as before. If ft; C ft - up then the trellis C; is 
a subcode of the trellis C* up . As such, instead of using C; we 
can use C* up at the ith stage of decoding. ft* up can be chosen 
to have a simpler trellis structure as compared to that of ft*. 
This would reduce the complexity associated with closest coset 
decoding and hence reduce the decoding complexity associated 
with the ith stage of decoding. 

Multistage decoding leads to error propagation. To reduce 
the effect of error propagation, the first couple of decoding 
stages should be powerful. A special case of the decoding 
algorithm occurs for q = 2 and £2 — ^ 2 - If closest coset 
decoding at the first stage is carried out in a single-stage, 
then the overall decoding of the multidimensional code is 
also one-stage. If m is small, then one-stage closest coset 
decoding is feasible, however if m is large, multistage closest 
coset decoding could be adopted to reduce the decoding 
complexity. The overall decoding in the latter case would then 
be multistage. 

Decoding Complexity of the Proposed Decoding Algorithm: 
The complexity of the proposed schemes will be measured in 
terms of the number of computations required for the decoder 
to produce an estimate of each 2-D PSK signal. For 1 < i < 
let 7 i be the total encoder memory of the convolutional code 
used at the ith level in the proposed scheme. Consider the 
ith stage of decoding. Then, due to the Viterbi algorithm 
alone, the complexity is 2 ' Y * +fct additions and 2 7 i ( 2 *’ — 1 ) 
comparisons, per m x 2 -dimensions (since each branch has 
m MPSK signals). The branch metric calculation forms an 
additional complexity and depends upon the choice of the 
inner codes. Let us call this complexity Bc x . Hence the total 
complexity per mx 2 dimensions is: 1 ) E ^ =1 2 y ' +k ' additions, 
2) £? =1 2 7, (2 fci - 1) comparisons, and 3) Ef =i 1 3c, • Dividing 
this total complexity by m would give us the number of 
computations required per two dimensions (i.e., the number 
of computations required to decode a single MPSK point). 

VI. Design Rules For Good Codes 

The performance of codes designed using the proposed 
technique depends upon various factors. If all the design 
considerations are followed strictly, the codes usually would 
achieve good performance and in some cases, with reduced 
decoding complexity. Some of the most important design 


considerations are: 1 ) the number of levels q , in the multi- 
level concatenation should be kept as low as possible. The 
advantages of this are twofold. First, reducing the number of 
encoding levels, would reduce the number of decoding stages 
and in most cases reduce the decoding complexity. Second, 
reducing the number of decoding levels also decreases the 
amount of error propagation which occurs as a result of the 
multistage decoding. To reduce the error propagation due to 
multistage decoding, the first few levels should be chosen 
extremely powerful, so that the amount of error propagation is 
decreased. This however leads to higher decoding complexity 
for the first few levels. 2) The number of dimensions, i.e., 
m x 2, should be kept as low as possible. As m increases 
the number of nearest neighbors associated with the code 
also start increasing, which limits the performance of the 
code. On the other hand, increasing m usually helps in 
decreasing the normalized decoding complexity associated 
with the code. 3) Theorem 5 gives us the minimum squared 
Euclidean distance of the overall multidimensional TCM code. 
For a given minimum squared Euclidean distance of the 
TCM code, de- free of the convolutional codes chosen to 
form the multidimensional TCM code should be chosen to 
be as small as possible. Lower de-free would imply lower 
decoding complexity associated with the convolutional code 
decoding. The above also holds for Theorem 6 . 4) The branch 
computation complexity Bq at the ith stage of decoding 
depends upon A*. If A* is chosen to have a simple trellis 
structure, the corresponding branch computation complexity 
will be minimal. If on the other hand, the trellis for A* 
is sufficiently complex, techniques described in Section V 
can be used to reduce the computation complexity. These 
techniques however, usually lead to degraded performance. 
5) Construction of codes with good phase invariance, places 
restrictions on codes as per Theorem 7 and hence in most cases 
this would limit either the performance and/or the achievable 
spectral efficiency. 

Most design considerations mentioned above lead to con- 
flicting requirements. Hence, there is a tradeoff involved 
between performance, decoding complexity, spectral efficiency 
and phase invariance. 

Appendix 

Proof of Theorem 7: The proof follows very closely the 
derivation of the phase invariance conditions in [16]. For the 
code to be phase invariant by 180°/2*“\ any code sequence in 
the multidimensional code when rotated by 180°/2 £_/l should 
produce another code sequence. Let V be the transmitted code 
sequence. Let V TOt denote the code sequence V rotated by 
180°/2*“\ Recall from Section III, that the basic building 
block of the proposed multidimensional codes is Ao> hence 
any valid code sequence in the multidimensional code can 
be considered to be a sequence of points from Aq. Consider 
the jth time instant. Let V j — A (V 1 j * V 2 j * ■ * * V t j) 
be the transmitted sequence of m MPSK signals at the jth 
time instant, where, V^ j € Co,, for 1 < i < L Also, let 
Vj ot = * '*V™j) be the sequence of m MPSK 

signals for V TOt at the jth time instant, where, V\°j € Co,, 
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for 1 < i < i. Using results of [16], Vj ot can be written in 
the following form: 

Vj ot = A((V X j + n j) * (V 2j + r 7 j)* 

'-'tYtj + Kj)) (AJ) 

where V±j = 0 for 1 < i<h, Vfr j — 1 and j = 
Vfoj • v h+1 j- j for h < i < £ and 0 denotes the 

all-zero sequence of length m. Form the sequence V\ such 
that the j\h time instant of V f is 

= < A - 2 > 

Then, for the code to be phase invariant under rotations 
of 180°/2*“ h T V' should also be a valid code sequence. 
Sequences of signal points from form a valid code 

sequence. Hence, if V 1 - € then V is phase invariant 

under rotations of 180°/2* - \ i.e., if 1 € C q ~\,h and 

V hj V h+ij V i-iJ^ C o-^ for h<i<t (A.3) 

then V is phase invariant under phase rotations of 180°/2* -/l . 
Since the above should hold for any transmitted sequence V , 
the Theorem follows. 
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